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Abstract 

The Stieltjes coefficients 7fc(a) arise in the expansion of the Hurwitz zeta 
function C(s,a) about its single simple pole at s = 1 and are of fundamental 
and long-standing importance in analytic number theory and other disciplines. 
We present an array of exact results for the Stieltjes coefficients, including se- 
ries representations and summatory relations. Other integral representations 
provide the difference of Stieltjes coefficients at rational arguments. The pre- 
sentation serves to link a variety of topics in analysis and special function and 
special number theory, including logarithmic series, integrals, and the deriva- 
tives of the Hurwitz zeta and Dirichlet L-functions at special points. The results 
have a wide range of application, both theoretical and computational. 
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Introduction 

The Stieltjes (or generalized Euler) constants 7A;(a) appear as expansion coeffi- 
cients in tlie Laurent series about s = 1 for tlie Hurwitz zeta function ({s, a), one of 
tlie generalizations of the Riemann zeta function C('^)- Elsewhere [6], we developed 
new summatory relations amongst the values 7fc(a) as well as demonstrated one of the 
very recent conjectures put forward by Kreminski [18] on the relationship between 
7fc(a) and — 7fc(a + 1/2) as A; — cxd [5]. 

In this paper, we present an array of exact results for the Stieltjes constants. These 
include individual and summatory relations for paired differences of these coefficients 
for rational arguments. Our work provides a unification of several important topics 
of analysis and analytic number theory. These include certain logarithmic sums, inte- 
grals of analytic number theory, special functions, the derivatives of the Hurwitz zeta 
and Dirichlet L-functions at special points, and differences of the Stieltjes constants. 
As the corresponding logarithmic sums are slowly converging, our results provide use- 
ful complements and alternatives to numerical computation. In addition, our analytic 
results and accompanying representations of 7fc(a) provide a basis for inequalities and 
monotonicity results for the Stieltjes constants. 

We may stress the fundamental and long-standing nature of the Stieltjes coeffi- 
cients. They arise in the expansion of the Hurwitz zeta function about its unique 
polar singularity. They can be used to write other important constants of analytic 
number theory, and they appear often in describing error terms and as a result of 
applying asymptotic analyses. As well, they can be expected to play a role in investi- 
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gations of the nonvanishing of L-functions and their derivatives along the hne s = 1 
and elsewhere in the critical strip. 

An expository paper of Vardi discusses the evaluation of certain logarithmic 
integrals (or their equivalents through change of variable) and describes the underly- 
ing connection with Dirichlet L-series. However, the presentation is illustrative and 
no connection with logarithmic series or the Stieltjes coefficients is mentioned. Much 
more recently, Medina and Moll [I9] have followed Vardi's approach and given a num- 
ber of examples for integrands containing a rational function. Additionally, Adamchik 
[2] considered differences of the first derivative of the Hurwitz zeta function at ratio- 
nal arguments and related them to logarithmic integrals. Below, we explicitly write 
one of his results in terms of a difference of values of 71(a). 

After exhibiting some basic relations for C{s,a) and 7fc(a), we provide both gen- 
eral relations to logarithmic sums and instances of logarithmic integrals in terms of 
the Hurwitz zeta function and iterated logarithmic integrals in terms of differences 
of Stieltjes constants. Motivated in part by such relations, we then proceed to other 
new analytic results, systematized as Propositions. We provide selected examples for 
differences of low-order Stieltjes coefficients. However, our methods apply to 7A;(a) 
for any order k for rational arguments a. We further obtain series representations 
of the Stieltjes coefficients. Among these results, we indicate how families of rapidly 
convergent summation representations may be obtained, that are applicable to com- 
putation. 



3 



Preliminary definitions and relations 

The Hurwitz zeta function, initially defined by 

oo 

can be analytically continued to the whole complex plane C — {1}, satisfying, for all 
s, the functional equation 



C(l - s,p/g) = ^E-s. - ^) as,r/g), (0.2) 

where p and q are integers and F is the Gamma function. Proof of the Hurwitz 
formula [25l EH] 

(27?)" , TTS ^ cos27ma . vrs sin27rna ^ ^ . 

i^Cil -s,a) = cos— + sm— ^ , Re s > 1, 0.3 

2r(s) 2 2 



has been simplified in References [3] and [30] . 

The defining relation for the Stieltjes constants in terms of a Laurent expansion 

is 

as,a)^^^t^-=lhM(s s^l. (0.4) 

S i AC. 

This equation reflects that ({s, a) has a simple pole at s = 1 with residue 1. By 
convention, 7^ represents 7a:(1) and thus explicitly we have [U [m [TTl [201 

1 00 / -1 \k 

C(.) = ^ + E - ^^1- (0-5) 

As noted in [29] and again in [30|, we have 70(a) = —il){a), where i/) = F'/F is the 
digamma function, so that 70 = 7, the Euler constant. A representation for the 
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Stieltjes constants is given by [3] 



7fc(a) = lim 

W— >oo 



^ ln^(j + a) _ \n^+\N + a) 



(0.6) 



(fc + l) 

As written, Eq. (0.6) is not well suited for numerical purposes. 

A very recently derived representation [8J (Proposition 3a) that is useful for com- 
putation and other purposes is given by 
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In'^+V 2 



7^(a) = — In'^a- + -Re / — ^ \ ^ 



dy, Re a > 0. (0.7) 



In the following, we let s{n,m) denote the Stirling numbers of the first kind and 
Ck{a) = 7fc(a) — (In'^ a)/a. With Bn{x) the BernouUi polynomials, their periodic 
extension is denoted Pn{x) = Bn{x — [x]). Two other integral representations of the 
Stieltjes constants of interest in this paper are [30] 



c„(a) = (-ir'd|:£<!i-+i."+i-*) 

k=0 



roo rr 

/ Pnix-a)^^dx, n>l. (0.8) 



kl 



and 



7n(a) 



E 

k=0 



ln"(A; + a) 
k + a 



ln"+i(m + a) ln"(m + a) p 

+ / Pi{x)fn{x)dx, 



n + 1 



2(m + a) 



< a < 1, 



0,1,2,..., 



(0.9) 



where fn{x) = ln"'(x + a)/{x + a). 
From Eq. (0.4) we have 



ds 



[as,a)-as,b)]\ 



s=l- 



(0.10) 
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Therefore, from Eq. (0.1) we obtain the relation to logarithmic sums 

/N /LN ^\ln\n + a) ln^(n + 6)" 

n=0 . 



(0.11) 



n + a n + b 

We may stress here a theme of this paper, namely that in Eqs. (0.10) and (0.11) 
polar singularities have been precisely eliminated. 

Various integrals may written in terms of differences of Hurwitz zeta functions. 
For example, we have 

a + r 



1 t'^-i In"-^ t 

and 

1 r-i In'-^ t 



dt 



^ Tis) 



1+tP 



dt^i-- Vis) 



C s, 



a+p 



Re s > 0, Re a > 0, 
(0.12) 

, Re s > 0, Re a > 0, Rep 
(0.13) 



Equation (0.12) is easily proved by writing the integrand factor {1 + t)~^ as geomet- 
ric series, interchanging summation and integration, and integrating termwise. The 
interchange is readily justified on the basis of uniform convergence. Equation (0.13) 
may be similarly obtained, or through the use of the change of variable u — V and 
the application of Eq. (0.12). 

Differentiation of Eq. (0.12) or (0.13) with respect to the parameters a, p, and/or s 
gives families of integrals. Specifically, differentiating Eq. (0.13) j times with respect 
to s gives 



L 



1 In"-^ t 



d 



1 



lnnntrft= - -- r(.) 



ds 



2p, 
6 



a + p 



Re s > 0, Re a > 0, Re p > 0. 



(0.14) 



In this way, we may relate such integrals to sums of differences of Stieltjes constants. 
As an example special case we have 



1 r^Mn Int , d 
l + tP ds 



— [7-i7r + ln(2p)] 



i( 2p) 

' a + p 



T{s) 



2p 



1 

2p 



[Tp 



2pJ 
a + p^ 



(0.15) 



A result of Adamchik [2] (Proposition 1) holds for integers p and q with p/q < 1. 
It states that 



p 



C 1, 



p 



71 1, 



p 



71 1, 1 



p 



TT cot 



[ln(27rg) +7] - 27r|]^lnr (^^j sin (^^^ . (0.16) 



This result was obtained on the basis of differentiating Rademacher's formula and is 
attributed to Almkvist and Meurman. 

Summatory relation for the difference of Stieltjes constants 

We have 

Proposition 1. We have the summatory relation 



-r[7n+i(a) - 7„+i 



n=0 



In 



m 

r(a) 



Re a > 0, Re 6 > 0. 



:i.ii 



Proof. We supply two different proofs. For the first, we apply the relation 
(pp. 87, 92) 

C'(0,a) = Inr(a) -lnV2^. (1.2 
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From Eq. (0.4) we have 
C{s,a) 



s - 



Ip + E^^7n+i(a)(.-ir. 



ra=0 



:i.3) 



Then putting s = and using Eq. (1.2) we have 



C'(0,a)-C'(0,6)=ln 



r(a) 

r(6) 



:i.4) 



and the Proposition follows. 

For the second proof we initially assume that a > and 6 > 0. We apply the 
representation (0.7) for 7fc(a) and have 



— 7yfc+i (a) = - In a + a — a In a — 1 H — Re 

K ! ^ 



fe=0 



(y/a - i)(a - i?/)ln(a - iy) 



a Jo (1 + ?/Va2)(e2'^f - 1) 



a+{--l-aj\na + 2lmJ^ j^-—-^dy, 



dy 



[1.5) 



where we have used elementary exponential sums, and Im ln(a — iy) = — tan ^{y/a). 
As j^ip{s)ds = ln[r(6)/r(a)], the Proposition follows from Binet's second expression 
for Inr(z) [23] (pp. 17, 91), 

\nT(z)= (z--)\nz-z + -\n{2n) + 2 r^-^^i-^^^dt, Re ;2 > 0. (1.6) 
V 2/ 2 JO e^'^* — 1 

By analytic continuation, the result (1.1) is extended to Re a > and Re 6 > 0. 

Remarks. The representation (0.9) or its equivalents may also be used to prove 
Proposition 1, in view of [9J (p. 107) 

r°o Pi{x) 



1 r°° Pi ( x) 

lnr(s + l) = (s + l/2)lns-s + -ln27r- / , ^ \ dx, Re s > 0. 
^ ^ ^ ' ' 2 Jo (x + s) 



[1.7) 



The representation (0.7) may be used to develop many other summations, includ- 
ine for \z\ < 1, 



^, ,.-ri(a) = -a'~Mna + (a"-l)---lna + 2Im -^-^——-^\n{a-ty)dy. 

(1.8) 



k=0 



Since C'(^l) — 1/12 — In A, where is A is Glaisher's constant, we find from Eq. 
(1.3) that 

OO 2'T' 1 

$:-7n+i = lnA--<0. (1.9) 

n=0 

The polygamma functions ip^^"^ are connected to the Hurwitz zeta function via 
7/>(")(x) = (-l)"+ir2!C(n + l,x) for inte gers n > 1. Therefore, we obtain from Eq. 
(0.4) for the trigamma function 

OO / -\\k 



k=0 



k\ 



-lk[.a) 



:i.io) 



and more generally 

^W(a) = (-l)"+in! 



1 - (-1)'= , 



n > 1. 



:i.ii) 



This equation may be taken as an infinite linear system. Its inversion would yield the 
Stieltjes coefficients in terms of polygammic constants. 

Relation to derivatives of Dirichlet L-functions 

We now introduce Dirichlet L-functions L±k{s) (e.g., [15], Ch. 16), that are known 
to be expressible as linear combinations of Hurwitz zeta functions. We let Xk be a 
real Dirichlet character modulo k, where the corresponding L function is written with 
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subscript ±k according to Xk{k — 1) = ±1. We have 

W^) = f:^ = ^Ex.Mc(^,y), Re.>l. (2.1) 

n=l ^ 1^ m=l ^ ^ 

This equation holds for at least Re s > 1. If Xfc is a nonprincipal character, as we 
typically assume in the following, then convergence obtains for Re s > 0. 

The L functions, extendable to the whole complex plane, satisfy the functional 
equations [31] 

and 

= i(27r)^fc-^+i/2 sin 1(1 - s)L^^{\ - s). (2.3) 

Owing to the relation 

r(i-s)r(.) = -^, (2.4) 

sm vrs 

these functional equations may also be written in the form 

L_,(l -s) = 2(27r)-^F-i/2 sin V{s)L^^{s), (2.5) 



and 



L+,(l -s) = 2(27r)-A;^-i/2 cos ) r(.)L+,(.). (2.6) 



Integral representations are known for these L-functions. In particular, we have 
Lemma 1 [31] . For a nonprincipal Dirichlet character we have 



1 /-oo U"'^ ' ^ 



L±j.is) = —-J^ ^—-^[J2xkim)e-'^-]du, Re s > 0. (2.7) 



\in=l 
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Remark. We may note that convergence obtains in the stated domain due to 
contribution in the numerator arising from cancellation due to the nonprincipal char- 
acter. 

Proof. For Re s > 1 and Re a > 0, we have the integral representation (e.g., [23] . 
p. 89) 



^s-l -(a-l)t 

Cis,a) = ^i _ ^ dt. (2.8) 



Tis) Jo 



We use Eq. (2.1) to write 



L±k(s) 



k- 



Jo 



oo j-s-l^-(m/k-l)t 



dt. 



(2.9) 



We then sum the finite geometric series and put u = t/k, giving the Lemma for Re 
s > 1. By analytic continuation for Xk nonprincipal, it also holds for Re s > 0. 

Key relations upon which we build are contained within the following. 
Proposition 2. For Xk a nonprincipal Dirichlet character we have (a) 



d_ 
ds 

k 



L±k{s] 



s=l 



E 

n=2 



Xk{n) Inn 



n 



k-'Exkim) [C (l,^)-lnA:c(l,^ 

k 



m=l 



1 — e 



^i\ — ]+\nk^/ji^j 
"^A du 



^ L Xk{m)e 



and (b) 



L±k{s] 



\m=l 



= E 

=1 n=2 



Xk{n) In^n 



n 



(2.10) 

(2.11) 
(2.12) 
(2.13) 

(2.14) 
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k ^ Xk{m) 



m=l 

k 



k ^ Yl Xk{m) 



m=l 



°° [ln2'U + 27lnM + 72 -((2)] 



(2.15) 
(2.16) 
(2.17) 



Kni=l 



1 - e-^" 

Proof. For part (a), Eqs. (2.10) and (2.11) follow from Eq. (2.1). For Eq. (2.12), 
we differentiate Eq. (2.1) and use the important relation for Xk nonprincipal 



ExfcW = i:x.(r) = 0. 



(2.18) 



r=l r=l 

Evaluation of the result at s = 1 gives Eq. (2.12). For Eq. (2.13) we have from 
Lemma 1 



d 



L±k{s) 



1 /-^ ii*~^[lnu - ■0(s) 



r(s)Jo 1 



3 — feu 



E Xk{m)e-^A du, (2.19) 

.m=l / 



where we used r'(s) = T{s)ip{s). Evaluating at s = 1, with = —7, gives Eq. 
(2.13). 

Part (b) follows very similar steps. In obtaining Eq. (2.17), we use ^'(1) = C(2) = 
7r^/6, where ip' is the trigamma function. 

Remark. A key feature of Proposition 2 connecting the values -£/±fc(l) with 
differences of Stieltjes constants is the nullification of polar singularities. An example 
of various relations of this Proposition, Eq. (0.11) at ^ = 1, and Eq. (0.15) at p = 1 
is given by 



a + 1 



d_ 
ds 



ry 



(2.20a) 



12 



= ln2 



When a — 1/2, the Dirichlet L-function appearing on the right side of Eq. (2.20a) is 



^-4(5). 



Case of Xk{k — 1) = — 1 



We have 

Proposition 3. Suppose that Xk is a nonprincipal character and that Xk{k — 1) — —1. 
Then 

Yl XfeMTi -r = ^^-fe(l) - In ^ Xk{m)ip I — 

m=l ^ ^ m=l ^ 



= k 



°° (\nu + 7) 
1 — e~*^" 



\m=l / m=l 



(ln27r + 7) E mx.(m) - TrA;^^ in JJ r^"^""^ (?) "^^^ E X.H^ 

m=l m=l V ft / j„=i \ 



(2.21) 

Proof. The first two equahties in Proposition 3 follow directly from Proposition 
2. For the third, we express L'_f^{l) in terms of known quantities of L'_^{0) from the 
functional equation (2.2) or (2.5). Given that C(O)C^) — 1/2 — a, we have from Eq. 
(2.1) that 

L±k{0) = E Xk{m) - ?) = E ^Xkim), 

m=l ^ '^m=l 

where we used property (2.18). Parallel to Eq. (2.11) we obtain 



(2.22) 



g-^L^k{s) 



k ^ E Xkim) 



s=0 



m=l 



C'(0.^)-..C(0.^) 



/m 



-Ink L±k{0) + J2 Xk{m)\nT [-) , 

m=l ^ 



(2.23) 
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where we used both Eqs. (1.2) and (2.18). From the functional equation (2.5) we find 

^1/2 ^1/2 

-L'_^(0) = (In A; - ln27r - 7)^-fc(l) + (2-24) 

TT TT 

This equation can be solved for L'_^(l), where Eq. (2.23) gives L'_j^{0) and by the 
functional equation (2.2) we have L_k{l) = (n /k^^'^)L^k{0)- The result is 

^'-.(1) = -^(ln2vr + 7) E -X.(m) - n^r-^-^ {j) , (2.25) 

yielding, by the first line of Eq. (2.21), the conclusion of the Proposition. 

We note that the digamma function, like the Gamma function itself, satisfies a 
number of identities that may be used to re-express Eq. (2.21). These include the 
reflection formula 

il){z) z) = -n coin z, (2.26) 
as well as the multiplication formula for integers m, 

ijj{mz) =\nm+ J2 ^ • (2-27) 

A;=0 \ "^z 

In addition, we recall that by means of Gauss's formula ([23], p. 19) the value 
ip{p/q) for any rational argument can be written as a finite combination of elementary 
function values. As concerns differences of higher order Stieltjes coefficients, similar 
identities may be written for the trigamma and higher order polygamma functions. 

Examples 

As examples of Proposition 3, we may write the following, using Eq. (2.26). 
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I ln27r + 7 - 31n 



r 



3yJ 



+ ln3^ 



(2.28) 



71 I 4 



-TT < In Stt + 7 — 2 In 



r 



(2.29) 



7r|2^(ln27r + 7) - Vein 



+ V61n6L (2.30) 



and 



7i(^;+7i^^ 



7H^ 



= V7n I (In 27r + 7) - In 



r 1 r f r | 



r 



r (f ) r (f ; 



7H^ 



ln7 , 



(2.31) 



[ji)-^' (ti)-^^ (n)-^^ (n)-^^ (n 



= TT 



^(ln27r + 7) - Vllln 



r {^) r ( r (a) r (a) r (a) 
r(^)r(^)r(^)r(A)r(if)J 
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In such equations, we could just as well use the duplication formula 



r{x) 



(2.32) 



(2.33) 



r(2x) ^ r{x+l/2) 

to re-express the Gamma function ratios. More generally, we may use the multipli- 
cation formula 



n-l 



T{nx) = (27r)(^-")/2^"^-i/2 ]J p 



X ■ 



k=0 



n 



(2.34) 
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Proposition 3 shows that the logarithmic sums (0.11) for rational values of a and h 
are essentially logarithmic constants. The integrals corresponding to the examples 
(2.28)-(2.32) are easily written from Eq. (2.21) and we omit the details. Comparison 
can be made to tabulated integrals that are expressible in terms of logarithmic ratios 
of Gamma function values [2[ [TTl [191 [26] . These include entries on pages 532, 571-573, 
and 580-581 of a standard table pT] . 

Regarding examples (2.28), (2.31), and (2.32) we may recall that for k an odd 
prime, there is exactly one nonprincipal Dirichlet character Xk modulo k. 

By [5] (Proposition 3) or as a special case of Proposition 5.1 of for integers 
g > 2 we have 

^"^ /r\ In^+^rt ^ fk\ 

E 7^ [-) = -7. + ^(-1)' + y (-l)^(l^^' ^)7^-r (2-35) 

Given Proposition 3, we may now combine various sums and differences of Stieltjes 
coefficients to find identities in terms of the values 7^ = 7^(1). As a very particular 
instance, we have 

Corollary 1. The values 7i(l/3) and 7i(2/3) may be separately written in terms of 
71- 

This statement follows from Proposition 3 at A; = 3 [Example (2.28)] together with 
Eq. (2.35) at g = 3 and = 1. 

In the case of 71(a), the relation (0.16) of Adamchik may also be invoked. 

Case of Xkik — 1) = +1 

We have 
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Proposition 4. Suppose that Xk is a nonprincipal character and that Xk{k — 1) — +1. 
Then 

m=l ^ ^ m=l ^ 



2(ln27r + 7) In fl r>^'=^"^^ (t) - E X^(HC" (o, ^) 



m=l ^ ' m=l 



m=l 



Proof. When Xfe(A; - 1) = 1, L+fe(0) = 0, E^i mxfe(m) = 0, and 



m=l 



Differentiating Eq. (2.1) we have 



L+k{0) = Yl Xk{m) 

m=l 



Differentiating the functional equation (2.3) we find 



m 
J 
(2.36) 



(2.37) 



(2.38) 



2(7 + ln(|^))LV,(0)-L'|,(0) 



(2.39) 



Then using Eqs. (2.35) and (2.36) we determine 



2(ln27r + 7)ln {[ T^^^^^ (^) - E Xk{m)C" (o,^ 

m=l ^ ^ m=l ^ 



Substituting into the first hne of Eq. (2.36), the Proposition is completed. 
Examples. We have from Proposition 4 at A; = 5 using Eq. (2.4), 



(2.40) 
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x/5 {2(7 + ln27r) In \{l + Vh) - (" (o, J) + C" (o, ^) + C" (o, ^) - C" (o, ^ 

+21n5coth-^ v^} . (2.41 



We have from Proposition 4 at A; = 10 using Eq. (2.4), 

3 \ / 7 \ / 9 



^'^loj'^' [lo]-^' lioJ+^Hio 



10 {2(7 + m 2.) m 1(3 + v/5) - c" (o, 1) + C" (o, ^) + c" (o, ^) - C" (o, ^ 

+3v/2 In 10 coth"^ Vs} . (2.42) 
From the well known Hermite formula for C{s,a) (e.g., [23], p. 91) 



a 



l-s 



oo 



sin (s tan ^y/a) dy 



we obtain 

("(0, a) = 1 In^ a + 2a In a - a In^ a - 2a - 2 tan"^ 

Evaluation of the integral for rational values of a, or combinations of rational values 
of a, would be of interest in regard to Proposition 4. It appears that a contour 
integral evaluation may be possible, with the integrand having simple poles along the 
imaginary axis at y = ji and residues (i/27r) tanh~^(j/a) ln(a^ — j^) there. However, 
this is likely to give an infinite series representation of C"(0) whereas another closed 
form is desirable. 

An explicit summation expression for the Stieltjes constants 

In this section we develop an explicit formula for 7n(a), based upon the represen- 
tation (0.9). What is required is an evaluation of the logarithmic integral there. We 



have 



Proposition 5. Let r(s,t) be the incomplete Gamma function. Then for < a < 1, 
m, n = 0, 1, 2, . . ., we have 

ln"(A; + a) ln"+^(m + a) ln'^(m + a) 



7n(a) = 



k=0 



k + a 



n+1 



2(m + a) 



+ ln"(j+a + l)-ln"(j + a) 



1 



ln"+^(j + a + l)-ln"+^(j + a 



{n + l) 

-(a + J + 1/2) [r[n, ln(j + a)] - V[n, ln(j + a + 1)] - T[n + 1, ln(j + a)] + V[n + 1, ln(j + a + ] 

(3.1) 

Proof. We may proceed as follows: 
/ Pi{x)fUx)dx=J2 Piix)jUx)dx=J2 (x-j-l/2)/;(x)rfx. (3.2) 



By using the expression for /„, we obtain 



/ p,ix)ax)dx=Y: /. 

n-1 



--(a + j + 1/2)- 



X 



x" 



{n In" ^ X — In" dx 



E -(« + i + V2)/ ^ . -dx- 



3+a X^ {n + 1) 

+ ln"(j + a+l)-ln"(j + a)}. 



ln"+'(j + a + 1) -ln"+^ (j + a) 
(3.3) 



The remaining integrals may be performed as [TT] 

fi+a+i In" X 



dx = V[n + l, ln(j + a)]-V[n + 1, ln(j + a + 1)]. 



(3.4) 



/j+a X^ 

The insertion of Eqs. (3.3) and (3.4) into Eq. (0.9) provides the Proposition. 

Remarks. Semi-infinite integrals over Pi as we have just performed are of much 
interest in connection with applications of Euler-Maclaurin summation. 
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In conjunction with Eq. (3.1) we may note that for n a nonnegative integer we 
have [nj (p. 941) 



r(n + l,x) =ri!e-^ 



n 



m=0 



ml 



(3.5) 



We may obtain other exphcit summation representations of the Stieltjes constants 
by working with Eq. (0.8). Here we show the approach for the special case of 71. 
However, by using the defining relation Bn{x) = Z]fc=o (fe)-^fc^""^' where Bj are 
Bernoulli numbers, thereby introducing another sum, this can be carried out generally. 

As an example, we have 

In*' a; 



k=0 



k\ 



Pi{x) 



-dx. 



(3.6) 



where s(2, 2) = 1, s(2, 1) = -1, and s(2, 0)/2 = 0. The integral of Eq. (3.5) is given 
by 



J I X 7 = 1 "'J 



Z-(i + l/2)i' 



X 



x^ 



In^ xdx 



E 



\n'^\3 + l)-\n'+'3 -(j + 1/2) 



k+l 



j+1 In'' X 



dx > . 



x^ 



(3.7) 



Denoting the integral = n~^^{\n'^ x)/x'^dx, we find by integration by parts the 



recursion relation 

I,{k) = kI^{k-l) + {j + l/2) 



J + 1 J 



(3.8) 



The solution of the homogeneous part of this recursion relation is of course just 
k\lj{0). Another way to evaluate the integrals Ij{k) is in terms of the incomplete 
Gamma function, as in Eq. (3.4). As Eq. (3.8) indicates, when effectively a — >^ in 
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Eq. (3.4), Ij{k) can be expressed explicitly for any fixed power k. The combination of 
Eqs. (3.6), (3.7), and (3.8) gives an explicit summation representation of the constant 
7i ~ -0.0728158454836767. 

New expressions for 71(a) and 7jk(a) 

We have 

Proposition 6. Let Hn — Z]fc=i 1/A; be the nth harmonic number. Then we have for 
Re a > 



71(a) = --ln\a+l)+Y: 7^— 4T[(C(A:+l,a)-a-('=+i))//fc+C'(^+l,a)]+lnaln 
^ k=i + ^) 



Proof. By Eq. (0.6) Bit k — 1 we have 



a 

(4.1) 



7i(a) = lim 



ln(j + a) ln^(iV + a) 



lim 



^ Hj + a) _ _ 

j + a Ji-a x + a 

ln(j + a) ln(a; + a) /""^ ln(a; + a) 

i=o 



\n{x + a) 



dx 



ln(j + a) /•! ln(a; + a) r 
"f, j + a Ji-a X + a Ji 



X + a 



dx 



= _1 ,„V + 1) + 5^ + E ^^^^ - E r i^^^rf- 



^_ ln2 a + 1 + + ^ / 

z a ~[ Jo 

We now apply the generating function for harmonic numbers 

ln(l + z 



ln(j + a) ln(x + j + a) 



j + a 



X + j + a 



dx. 



(4.2) 



n=l 



1 + ^ ' 



-^1 < 1, 



(4.3) 
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to expand the integrand in this equation. We have 

Iny \n{x + y) ln{l+x/y) 
y X + y X + y 

lii(l + x/y) \ny 



\ny 



1 



x + y 



y 



i^ + x/y) 



x + y y^ 
- 1 



lln(l + x/y) Iny 
y i + x/y 



E(-i)' 



,X" 



k=l 



y 



k 



k=i y 

We obtain for the integral in Eq. (4.2) 



k+l 



X . 



(4.4) 



'dx = J2 



(-1)^= [Hk-Hj + a)] 



•lo ti U + fi + 1) U + aY+' 



(4.5) 



Summing this expression over j = 1 to oo gives 



OO ^ -j^^fc 



S + 1) 



.=1 + 1) ' 

[(C(A; + 1, a) - a-('=+^))//fc + C'(A; + 1, a) + a-('=+^) In a], 



(4.6) 



where we used the relation C(s, a+ 1) = C(s, a) — a"'' and its derivative. Inserting Eq. 
(4.6) into Eq. (4.2) gives the Proposition. 

We mention a second method for proving Proposition 6, by using the integral rep- 
resentation (2.8) for C(s, a). This method can also provide an integral representation 
of 7i(a). In particular, since 

X j-oo ^feg-(a-i)t 



C'(^ + l,a) = -/^ -^Intdt, 



(4.7) 
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we have the term of Eq. (4.1 

oo f ^-^\k 



S (k + 1) 



Cik + l,a) 



g-(a-l)i 

e* — 1 



t 



- 1 



t 



\nt dt. 



(4.8) 



This integral may be evaluated per logarithmic differentiation of 



1 

Jo 



e* - 1 



t t 



rft = r(/5)[a^-/5C(/3 + l,a)]. 



(4.9) 



Then one may apply the operator [d/d[3)p=Q, use the functional equation of the 
Hurwitz zeta function, and the relations C(0, a) = 1/2 — a and (1.4). 
For the other summation term in Eq. (4.1), using (2.8), we have 



i:^C(*+l.a)ft 



'{a+l)t 



l-e-*t 



7(1 + e*) + 2e*/2 cosh (0 In t 



du 



i/=-l/2 



(4.10) 



where is the modified Bessel function of the first kind (e.g., pjj, pp. 958, 961). 

Note 1. By integrating the generating function relation (4.3) we obtain the term 
of Proposition 6 



oo (T\k 



a + 1 



(4.11) 



Corollary 2. We have 



71 = -2 2 + E (^[(C(A: + 1) - l)H, + Cik + 1)] 

oo (_l\k 



(4.12) 
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Corollary 3. Proposition 6 permits the recovery of relation (2.35) at A; = 1, 



E7i(^) =(g-l)7i-gQ + 7)lng. (4.13) 



In verifying this statement we use 

Ec(a: + 1,^^ =(g^+i-l)C(A: + l). (4.14) 

Rapidly converging expressions for 7fc(a) for a complex domain of a may be ob- 
tained by using very recently constructed series representations of the Hurwitz zeta 
function [7]. As an instance of this we find 

Proposition 7. Let the generalized harmonic number H^^ = J2k=i ^/k'^ = ^"'(1) ~ 
ilj'{n + 1) = 7r^/6 — tp'{n + 1). Then we have for Re a > 1/2 

1 / 1\ °° 1 
71 (a) = --ln'[a--j + ^^^^^^ [H^kCi^k + 1, a) + C'(2A; + 1, a)] , (4.15) 



and 

-72(a) = ^ In^ (« - ^) + E ^k^2k + l) [^^'^ " ^2?)^(2A; + 1, a) 

+2H2kC{^k + 1, a) + C{2k + 1, a)] . (4.16) 
More generally, for k > 1 and n > 1 let 

r.(fc,m)^(-l)— + 0<m<n, (4.17) 

where {z)a = T{z + a)/r{z) is the Pochhammer symbol and s{n,m) are the Stirling 
numbers of the first kind. Then we have 

ln{a) = ln"+i (a--) -(-1)" ^ E ^ m)C(™)(2A;+l, a). 



n+1 \ 2J ' ' ^,A\2k + l)-,\m, 



(4.1^ 
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Proof. We use [7] (Corollary 3) for Re a > 1/2, 



(a-l/2)^-% (a -1/2)1- 

C (s, a) = ^ ln(a - 1/2) - 



s - 1 



s-1] 



-E 



[shk 



{[^(s + 2A;) - ^(s)]C(s + 2A;, a) + C'(s + 2A;, a)} 



,-4'=(2fc + l)! 

and apply the definition (0.4). For Eq. (4.16) we use the property 

^(s)pfc = {s)pk[ip{s + pk) - ^/'(s)], 



(4.19) 



(4.20) 



so that 



ds 



{s)2k = {'2k)\H2k. 



More generally, we use the expansion 



T \ (_-\ )j . / 1 • 



(4.21) 



(4.22) 



2 / ^7 V 2 

and apply the product rule to find Eq. (4.18). 

New expression for 7n(a) 

Based upon Proposition 6, we expect that the corresponding summation expres- 
sion for 7fc(a) will contain a series of terms to order k. We have 
Proposition 8. Let s{n, m) denote the Stirling numbers of the first kind. Then we 
have for Re a > (i) 



1. 3,_^,,ln2_a_g 1 



72(a) = --ln^(a+l)4 



lyC'ik + 1, a + 1) - -s{k + 1, 2)C{k + 1, a + 1) 



k\ 



+-s{k + l,3)C{k + l,a + l) 
kl 



(5.1) 
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Our expectation is fully borne out, as our proof of part (i) extends to general n (ii): 

7n(a) = 



n + 



-ln"+^(a + l) + 

1 a 



OO -1 



(-l)^C^")(fc + l,a + l)--^ 



ni {-ly 



(5.2) 



Proof. We will use a generating function for the Stirling numbers s{j, k), 
1 



X 



ln'"(l + x) = m! ^ s(n + l,m)- 



X" 



\x\ < 1. 



(5.3) 



Proceeding similarly to the beginning of the proof of Proposition 6, we have generally 



7fc(«) 



k + 



-In^ a+1 + + E / 

1 a Jo 



ln''(j + a) \n''{x + j + a) 
j + a 



X + j + a 



dx. (5.4) 



Now we apply the expansion (5.3) to write for part (i) 



\n^{x + y) 



y 



x + y 



-E 

k=l 



y 



fc+i 



-l)'^ln^y + frs(fc + 1,2) In I/ + f-s(fc + 1,3) 



(5.5) 



We then substitute this equation into Eq. (5.4) at A; = 2 and perform the integration. 
We then perform the summation over j, using 



:-l)^C^^)(fc + l,a+l), 



and part (i) follows. Similarly, for part (ii), in place of (5.5) we use 



(5.6) 



ln"y ln"(x + y) 
y x + y 



E' 

k=l 



,k+l 



n-1 



_1)^ ]^ny^ J2 7z--—^s{k + 1, J + 2) W^-^-' y 



(5.7) 
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We substitute this equation into Eq. (5.4), perform the integration, perform the 
summation over j using (5.6), giving Eq. (5.2). The Proposition is complete. 

Remarks. Proposition 6 corresponds to the well known special case of the har- 
monic numbers wherein s{n + 1, 2) = (— l)"'^^n!if„. The connection of Stirling num- 
bers with sums of generalized harmonic numbers H^^ is well known. For example, 
we have s{n + 1, 3) = (-l)"n![iJ2 _ hP]/2. 

As for 7i(a), an integral representation for 7n(a) niay be developed using Eq. 
(2.8). We omit such consideration. 

If desired. Proposition 8 may be used to write an expression for the differences 
7„(a) — 7„(6). We have been informed that R. Smith has also obtained Eq. (5.2), and 
studied its rate of convergence [22] . 

Integral expression for 1 — 7 — 71, and more 

Let {x} = X — [x] denote the fractional part of x. It is known that 

/i = ^'|^|dx = l-7, (6.1) 

a result that is reproved and discussed in the Appendix. We show here 
Proposition 9. We have 

h = /' /' | — I dxdy = r r ^dydx = 1-7-71. (6.2) 
Jo Jo [xy ) Ji Jx xy^ 

Although this result is supposedly proved in Ref. [2T], we have not been able to 
obtain it, and so present our own proofs. These methods of proof may themselves be 
of independent interest. 
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Proof 1. We first establish 
Lemma 2. We have 



/-OO fy\ 

/ -^dy = Hi^] - 7-lna; + l- 
Jx y 



We have 



{y}. 



7^ y^ J[x] y2 



y' 



In 



-J 



[^] y^ 



(h) ^ 1^' (m 



The sum here is given by (cf. the Appendix) 



E 



In 



J + 1 



J + 1 



E 



In 



J + 1 



J + 1 



N-i 

-E 



In 



1-7-E 

i=2 



Inj - ln(j - 1) - - 
3 



= 1 - 7 - ln[a;] + H^^^] - 1 = if^^;] - ln[a;] - 7. 

Inserting this expression into Eq. (6.4) we obtain the Lemma. 
Using the Lemma, we have 



f°° f TT 1 . N \ dx 
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lim 

M— »oo 



1 ^ /-i+i / 

(l-7)lnM--ln2M + ^ / (h,- 
^ ,=1 Jj V 



M^oo 

Now the term 



= lim ^ (l-7)lnM- -In^M + y 



M 



1 



J + l \ 

. 3 ) J+1 



(6.6) 



fe=i ^ j=fe 



^-[ln(M+ 1) -InA;] 
fe=i ^ 



^ InA- 
^i^MMM+l)-^^. 

fc=l 



(6.7) 



Then by Eq. (6.6) we have 



r 1 ^ InA- 1 

h= hm Ul-7)lnM--ln'M + i/Mln(M+l)- V — -i/M+i + U. 
M^oo [ 2 A; J 

(6.8) 

We recall the asymptotic form Hm = InM + 7 + 0(1/M) as M — > 00 and appeal to 
Eq. (0.6) at a = A; = 1: 



72= lim |lnM-7lnM- -ln^M+ [7 + lnM + O 



^ A; 



7 + lnM + 0(^ 



M 

+ 1 



InM + of-J- 



lim i;Tln^M-V — -7 + 1 + —— 



1-7-71- 



(6.9) 
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The alternative form of I2 in Eq. (6.2) results from the change of variable {u, v) 
{l/xy,x), with inverse transformation {x,y) = {v, 1/uv), and Jacobian 



d{x,y) 




1 


1 


d{u, v) 




]_ ]_ 





Then we obtain 



1 roo 



{u] 



dudv 



00 roo 



{4 



dzdt. 



(6.10) 



(6.11) 



/O Jl/v U^V Jl Jt Z^t 

For the second equality we have used the simple transformation {v, u) = {l/t,z), with 
corresponding Jacobian 

' (6.12) 



d{u, v) 



1 



We may obtain a second shorter proof of Proposition 9 by using a known result, 
namely the integral representation [16] 

roo ]^ 



jlnH dPl{t)= ——[k-\Yit)P,{t)dt-5ko/2, 



Ik 



(6.13) 



where 5jk is the Kronecker symbol. We interchange the order of integration in the 
second form of I2 on the right side of Eq. (6.2), giving 

r-y {y} 



1 J I xy^ 



;dxdy 
If 



1 y 

Inv dv 



— Iny dy 



(6.14) 
(6.15) 



Therefore, from Eq. (6.14) and Eq. (6.11) at A; = 1 we obtain 

[Piiw) + l/2] 



-7 - 1 + 



1 w 

°° Pi{w) 
1 w'^ 



dw -\ — / 
2 J I 



\nw dw 

Inw 



dw 



w^ 
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= -7i - 7 + ^ + ^ = -7i - 7 + 1- (6.16) 

Remcirks. Since < {x} < 1 for x > 0, we easily have from Eqs. (A.l) and 
(6.11) or (6.14) the inequahties < Ji < 1 and < /2 < 1. 

We see from Eq. (6.13) that higher order Stieltjes constants may be obtained from 
higher dimensional integrals. For instance, there is a 72 contribution in the integral 

POO rz rz /rl 1 roo poo pz S z\ 

/ / / ^—dxdydz= / / / \^dxdzdy, (6.17) 
Ji Ji Ji xy Ji Jy Ji z^xy 

and such integrals may be further rewritten with changes of variables. 
Indeed, we introduce the integral 



h 

and have 



r r — \ dxdydz = r C t ^^dwdvdu, (6.18) 
^0 Jo Jo [xyz ) Ji Ji/u Ji/v w'vw 



73 = 1-7-71-^72. (6.19) 



This is a special case of the following. 
Proposition 10. Put 



In = [ / ■■■/ I \ dxidx2 ■ ■ ■ dxn, n>l. (6.20) 

Jo Jo Jo I XiXo ■ ■ - Xr, ) 



Then we have 



Corollary 4. We have 



4 = 1-E|. (6.21) 
j=o J- 



hm 7„ = -C(0) = I. (6.22) 
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Proof. We use the change of variable (xi, X2, . . . , Xn) = {un, Wn-i, ■ ■ ■ ,U2, 1/miM2 ■ ■ ■ Un), 
with Jacobian 



J 



■■■ 1 
■■■ 1 







1 ; 

1 ■■■ 



1 











U{U2 ■■■Un 



(6.23) 



where e„ = ±1. Specifically, if n is of the form 4m or 4m + 1, e„ = — 1, and e„ = 1 



otherwise. Then 



{^l} ni=n C^^i 

We then multiply integrate to find 



1 



[n — 1)! Ji u\ 



^In^-^Mi dm. 



(6.24) 



(6.25) 



We now note from Eq. (6.13) 



n-l 

E 

A;=0 



k\ 



1 



(n-l)! 7i 



r(n,lnt) dPi{t) 



Pi(t)cit+^, (6.26) 
[n — ly. Ji t'^ 2 

where we applied [11] (p. 941), integrated by parts, and used r(?T,,0) = (n — 1)! for 
n > 1 an integer. We then have the Proposition, as 



n-l 

E 

fc=0 



Jk 



1 



, , In"-^ t , 
{t}—:^dt+l. 



Remarks. We have found that the subjects of Propositions 9 and 10 have re- 
cently been of interest elsewhere [10]. It would be interesting to have a probabilistic 
argument for Corollary 4. 



32 



Summary 

We have obtained new explicit analytic results for the Stieltjes coefficients includ- 
ing series representations and summatory relations. Other integral representations 
based upon the properties of Dirichlet L-functions provide the difference of Stielt- 
jes coefficients at rational arguments, and these give inequalities. Our results have 
implications for other coefficients of analytic number theory and other fundamental 
mathematical constants. 
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Appendix: Integral expression for 1 — 7 



We here show that 
Proposition Al. We have 



h = N-]dx= [ 

Jo Ix) Jl 



1 X 

Proof. We recall that {x} = Pi{x) + 1/2 = x — [x], giving 

{x} 

■ux = 



(Al) 



(- - -^]dx 

^Jj \x x^ 



E 

N 



In 



J 



1 



lim 



j=l L 



ln(j + 1) - Inj 



J + 1 



lim [InN - Hn + 1] = 1-7, 



(A.2) 



where we used the telescoping nature of the sum in the next-to-last line. 

As a variation on this proof, we may note that a sum above is a case of the 
summation [12] (44.9.1, p. 290) 

k + x\ : 



E 

k=l 



In 



We then obtain 



k 



E 



In 



J + 1 



—jx + In r(a; + 1). 



(A.S) 



E 



In 



\ J J J+1 
'j + l\ 11 



J 



J J J + 1 



-7 + lnr(2) + 5: 



1 



1 



,=i\j J + 1, 
34 



1-7. 



(A.4) 



Similarly, this sum may be found at x = in [12] (44.9.4, p. 290) or at x = 1 in [12] 
(44.9.5, p. 290). 
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